On the superfluid properties of the inner crust of neutron stars 
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We investigated the superfluid properties of the inner crust of neutron stars, solving the Hartree- 
Fock-Bogoliubov equations in spherical Wigner-Seitz cells. Using realistic two-body interactions in 
the pairing channel, we studied in detail the Cooper-pair and the pairing-field spatial properties, 
together with the effect of the proton clusters on the neutron pairing gap. Calculations with effective 
pairing interactions are also presented, showing significant discrepancies with the results obtained 
with realistic pairing forces. At variance with recent studies on finite nuclei, the neutron coherence 
length is found to depend on the strength of the pairing interaction, even inside the nucleus. We 
also show that the Wigner-Seitz approximation breaks down in the innermost regions of the inner 
crust, already at baryonic densities pt > 8 • 10 +13 g ■ cm~ 3 . 

PACS numbers: 26.60.Gj, 21.60.Jz, 21.30.-x 



I. INTRODUCTION 

The inner crust of neutron stars [l[ offers a unique op- 
portunity to test nuclear-structure models. Its extremely 
neutron-rich environment represents a very strong chal- 
lenge to theories that have been developed for finite nu- 
clei. In particular, the Energy-Density Functional (EDF) 
method is able to describe to good accuracy properties 
of medium-mass to heavy nuclei 0, HJ and it has been 
applied to the inner crust of neutron stars since the sem- 
inal work of Negele and Vautherin [|| . They divided the 
inner crust of neutron stars into independent spherical 
Wigner-Seitz (WS) cells each of them representing a 
inner-crust region of a given density. 

To deal with a large number of nucleons, the EDF 
method relies on effective interactions that are fitted also 
to stable-nuclei experimental data. Recently, indepen- 
dent studies have been carried out by different groups 
to improve the connection of EDF theories to basic 
nuclear forces. Particular attention has been paid to the 
pairing correlations that are responsible for the superfluid 
properties of the nucleus. These studies use, in the pair- 
ing channel, phase-shift-equivalent interactions (so-called 
realistic interactions) evolved to low momentum through 
Renormalization-Group (RG) techniques [Iol - [l3l ]. 

The RG evolution is helpful in many respects. Not only 
it softens the original hard-core interactions, making the 
calculations feasible, but also it helps increase the EDF 
calculation reliability. This is because the hard-core po- 
tentials connect high-energy states, whose effective mass 



is not well described by the phenomenological EDF func- 
tionals, leading to reliability issues [l4j ]. Low-momentum 
interactions do not probe states too much high in energy. 

Along these lines, in this work we study the super- 
fluid properties of the inner crust of neutron stars, adopt- 
ing the WS approximation and solving the Hartrcc-Fock- 
Bogoliubov (HFB) equations with low-momentum real- 
istic interactions (V\ ow k) in the pairing channel. The 
pairing interaction matrix elements arc computed at 
first order only, leaving higher-order correlations for fu- 
ture investigations. For comparison, we also perform 
calculations with effective pairing interactions, namely 
the Gogny Dl interaction I15| and a Density-Dependent 
Delta Interaction (DDDI) [lj|, which have both been 
used in the past as pairing interactions in WS calcula- 
tions [T^-tHl- The superfluid properties obtained with 
these effective pairing forces turn out to differ substan- 
tially from those obtained with realistic pairing poten- 
tials (see Sects. [iTTBl and IjTTD]) . 

A comparison with the full band theory [lH, [23[ has 
shown that the WS approximation can reproduce well 
ground-state properties of the outermost regions of the 
inner crust. Its validity in the regions closer to the star 
core, where the clusters nearly touch each other, is still 
under debate. In this work we find that the WS ap- 
proximation starts to break down at baryonic densities 
Pb ~ 8 • 10 13 g ■ cm" 3 , where the protons leak out of the 
center of the cell and boundary-condition effects start 
kicking in. See Sect. Ill Al for a detailed discussion (see 



also Ref. j2j) 
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We study 11 different baryonic-density regions of the 
inner crust, from pb ~ 5 ■ 10 11 g • cm~ 3 ps 0.0018po to 
Pb ~ 10 14 g • cm~ 3 w 0.35po 5 with the saturation density 
pa = 2.8 • 10 14 g • enr 3 = 0.16 fm~ 3 . The correspond- 
ing WS-cell properties are shown in Table (|I| and have 
been taken from previous Hartree-Fock (HF) energy- 
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minimization calculations [4J. Recently. Hartree-Fock- 
Bardeen-Cooper-Schrieffer (HFBCS) |o|, [H[ and HFB 
fl9l ] minimization procedures have been carried out, ob- 
taining (Rws, Z) configurations that differ substantially 
from those in Ref. The actual cluster configurations 
in the inner crust still represent an open question and the 
energy-minimization calculations are very sensitive to the 
functionals used. However, as we show in this work, the 
supcrfluid properties of the system turn out to be rather 
independent of the (Rws, Z) configurations adopted for 
a given density region. 

In Sect. HJ we present the details of the calculations, 
followed by the results in Sect. Mil Sect. MI Al deals with 
the density profiles and in particular with the instabil- 
ities of the proton density for regions close to the star 
core. The effect of the proton clusters on the neutron 
pairing gap is discussed in Sect. IIIIB1 where results for 
Infinite Neutron Matter (INM) are compared to those in 
the inner crust. The spatial properties of the Cooper- 
pair wave function and of the pairing field are treated 
in Sects. MI CI and MI Dl Conclusions and outlook are 
presented in Sect. [IVJ 
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53.6 


4.67 


10 11 


0.12 
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6.69 


10 11 


0.15 
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250 Zr 


40 
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46.4 


1.00 


10 12 


0.19 
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320 Zr 


40 


280 


44.4 


1.47 


10 12 


0.23 
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500 & 


40 


460 


42.2 


2.66 


10 12 


0.31 
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950 Sn 


50 


900 


39.3 


6.24 


10 12 


0.43 


5 


110() Sn 


50 


1050 


35.7 


9.65 


10 12 


0.51 
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i350 Sn 


50 


1300 


33.0 


1.49 


10 13 


0.60 
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1800 Sn 


50 


1750 


27.6 


3.41 


10 13 


0.80 
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i500 Zr 


40 


1460 


19.6 


7.94 


10 13 


1.08 


1 


982 Ge 


32 


950 


14.4 


1.32 


10 14 


1.33 



TABLE I: The WS cells representing different density regions 
of the inner crust. The particle numbers Z,N, the WS-cell 
radii Rws and the baryonic density pi, have been taken from 
previous calculations [4j. kF, n is the Fermi momentum corre- 
sponding to the density of the outer neutron g £LS, ELS computed 
in this work. 



II. CALCULATION DETAILS 



1 TT 1 '1 
nlj 



A. Inner crust of neutron stars 

The self-consistent HFB equations [25[ are solved in 
each representative WS cell on a spherical mesh. The 
mesh step is 0.2 fm for the cells 11 through 8 and 0.1 fm 
for the higher-density cells (see Table (jH). The single- 
particle wave functions are expanded on a spherical 
Bessel basis with a momentum cutoff k max = 4 fin - . 
This corresponds to an HFB model-space energy cutoff 
of about H 2 k niax /2m « 320 McV. Our calculations are 
stable with respect to an increase of the model space and 
to a decrease of the mesh size. The WS-cell approxima- 
tion relies on the fact that the structure of the inner crust 
of neutron stars is recovered by a repetition in space of 
the WS cell. This requires the neutron density at the 
edge of the cell to be finite and to match that of the 
neighbor cells. This can be achieved by imposing the 
following Dirichlet-Neumann mixed boundary conditions 
[3| : (i) even-parity wave functions vanish at R = Rws', 
(ii) the first derivative of odd-parity wave funtions van- 
ishes at R = Rws- We call them Boundary Conditions 
Even (BCE), in contrast to the Boundary Conditions 
Odd (BCO) where the two parity states are treated in 
the opposite way. 

We use a Skyrme functional to build the single-particle 
Hamiltonian h and then we let the particles interact pair- 
wise in the pairing channel. The two-body matrix el- 
ements of the pairing interaction in the J = 0, T = 1 
channel enter the neutron-neutron and proton-proton 
gap equations, whose solutions provide the matrix ele- 
ments of the state-dependent gap matrix A. The latter, 
in turn, enters the HFB equations 



LS -nn'lj u n'lj /_^ l \ IL n'nlj fc F ,Q I v n> Ij ~ ^ilj V nlj 
n' n' 

(1) 

where ep^ q is the Fermi energy and q stands for neu- 
trons and protons. We used the standard notation nlj 
for the spherical single-particle states with radial quan- 
tum number n, orbital angular momentum I and total 
angular momentum j. U%j~j and are the Bogoliubov 
amplitudes for the z-th quasiparticle of energy Ef^ . 

When presenting the results for the HFB neutron 
pairing gaps, we show the Lowest-quasiparticle-energy 
Canonical State (LCS) pairing gaps [f|. The LCS gap 
is the diagonal matrix element of the gap matrix for the 
canonical state n a l a j a with the lowest canonical quasi- 
particle energy 



E n . 



laJa 



n a l a ja ~ £ F) 



A 2 



(2) 



where s na i a j a is the canonical single-particle energy. We 
dropped the isospin index q in Eq. [2] 

The Skyrme functional SLy4 has been used through- 
out this work, except for Figs. |4]an d[6l where a compar- 
ison with the functionals SkM*|H, [23 and MHF (Tj] is 
shown. 

We consider three different two-body pairing interac- 
tions: (i) a density-dependent contact interaction; (ii) the 
finite-range Gogny Dl interaction; (iii) low-momentum 
realistic interactions (Mowfc)- 

We restrict the Gogny Dl and the Vi ow & pairing inter- 
actions to the partial wave. Here we use the standard 
notation 2S+1 /,/ rol , with the Cooper-pair total spin 5, the 
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relative orbital angular momentum I and the relative to- 
tal angular momentum J le \ = l + S. Higher partial waves 
can also contribute to the superfluidity in finite nuclei, 
with P waves giving a « 15% quenching of the S-wave 
pairing gaps [jj . This contribution could be even smaller 
in the inner crust of neutron stars, where the states close 
to the Fermi surface are in the continuum and the center 
of mass of the Cooper pairs plays a less important role. 

We now give a detailed description of the pairing in- 
teractions that we used. 

(i) The two-body contact force DDDI between parti- 
cles at positions ri and r 2 reads 



v(ri,r 2 ) = V 



1-7/ 



IX) 



5{vi - r 2 ), 



(3) 



with V a = -430.0 MeV fm 3 , rj = 0.7, a = 0.45, 
Pa = 0.16 fm" 3 . We use a cutoff of 60 MeV 
on the quasiparticle energy. According to the 
literature [28j . this parametrization is such that it 
approximately reproduces the Gogny Dl pairing 
gaps in HFB calculations in INM. 

(ii) A separable pairing interaction that reproduces the 
1 So Gogny Dl pairing gap at the Fermi surface in 
INM (see Refs. [H, [30| for a detailed description) 

7j(ri,r 2 ,r' 1 ,r / 2 ) = 



(4) 



The operator |(1 — P a ) restricts the interaction to 
total spin S = 0. R = (ri + r 2 )/2 is the cen- 
ter of mass of the two interacting particles and 
r = ri — r 2 is their mutual distance. Strength and 
form factor are 7 = —738 MeV fm -3 and P(r) = 
l/(47ra 2 ) 3 / 2 exp(-r 2 /(4a 2 )), where a = 0.636. 



(iii) A rank-3 separable interaction 
form 



1 |3l|, |3J of the 



good approximation, with the neutron pairing gaps 
changing of at most 30 keV and 100 keV for INM 
and for the WS cells respectively, when the cut- 
off A ranges between 1.8 fin -1 and 4.0 fm" 1 . The 
G/3 (r) form factors are a product of a Gaussian and 
a Hermitc polynomial. 



Both neutrons and protons are found to be super- 
fluid in the WS cells of Table with the proton gaps 
comparable to the neutron ones. However, in this work 
we discuss only neutron superfluidity. We dropped the 
Coulomb term in the proton-proton pairing channel. In 
a few cells (i.e., 1350 Sn, 1800 Sn and 1500 Zr ), we checked 
that the inclusion of the Coulomb term in the Gogny Dl 
pairing interaction leads to a quenching of the proton 
pairing gaps between 20% and 30%. This is in agree- 
ment with recent studies on finite nuclei [33j . The neu- 
tron properties are not sensitively affected and discussion 
and conclusions are not changed by the inclusion of the 
Coulomb term, as neutron LCS gaps are affected at the 
level of 1 keV. 



B. Infinite neutron matter 

To study the effect of the proton clusters, the super- 
fluid properties of the inner crust are compared with 
those of the INM. For a given neutron density p n in INM, 
the HFB gap and number equations have to be solved si- 
multaneously 



A„(fc) 



d 3 k' 



v(k-k') 



An(fc') 

E n (k') 



Pn 



1 

2^ 



dk k 2 



1 



£n{k) /-In 



E n (k) 



(6) 



(7) 



p n is the neutron chemical potential and E n (k) = 
yj (e n (k) — p n ) 2 + A„(fc) 2 is the quasiparticle energy, 
while the single-particle energy e n (k) is given by the sum 
of the kinetic energy and the Hartree-Fock potential U^ F 



h 2 k 2 

2m* 



+ U n HF {k). 



(8) 



«(r 1) r 2j r' 1 ,r 2 ) = ^A^(r)G /3 (r')5(R-R / ) 

13=1 



(5) 



is used in this work to reproduce to high preci- 
sion the 1 So matrix elements of the low-momentum 
nucleon-nucleon interactions Vj w k obtained from 
the Argonne potential AV18 . The latter has been 
RG evolved to a low-momentum cutoff A using a 

6. The results shown in 



smooth regulator n 



exp 



this paper arc obtained with a separable force cor- 
responding to a A = 2.5 fm -1 low-momentum in- 
teraction. Our results are cutoff independent to a 



We use the Skyrme neutron effective mass m*. The 
number equation (cf. Eq. (|7"|)) provides the relation be- 
tween the density and the chemical potential p n - In the 
limit of weak coupling, where A n << EF.m the chemi- 
cal potential can be approximated by the Fermi energy 



n-ki 



-, with kpn = (37T p n ) ' . This approxima- 
tion somewhat holds already at fcjr n ea 0.2 fm -1 , and we 
are left with solving only the gap equation (cf. Eq. (|6|)). 
We check the validity of the above approximation a pos- 
teriori, by comparing the solution of Eq. (|6]) with that 
obtained in spherical-box calculations of homogeneous 
neutron matter, where both gap and number equations 
are solved simultaneously. The good agreement between 
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the two methods (see Sect. IIIIB|) supports our results. 
Since the neutron HF potential U^ F is constant and the 
single-particle energies are taken from the Fermi level, we 
also adopt the approximation e n (k) = f^r- 

III. RESULTS 

A. Density profiles and limits of the WS 
approximation 

HF and HFB neutron and proton densities for the 
inner-crust regions 11 through 3 (cf. Table (HJ) are shown 
in Fig. [1] In these regions, the HFB density profiles ob- 
tained with the three pairing interactions are almost on 
top of each other. The density of the outer neutrons 
gradually raises as one goes deeper and deeper into the 
inner crust, with the innermost WS cell 1800 Sn having an 
outer neutron density p n ~ 0.02 fm -3 . The proton clus- 
ters have a radial extension that ranges from 5 fm for the 
outermost cell (i.e., 180 Zr) to about 7.5 fm for 1800 Sn. The 
proton density above 10 fm is negligible for all regions 11 
through 3. 

At higher baryonic densities (pb > 0.25po)j the mean- 
field proton spatial distribution becomes unstable, as 
shown in Figs. [2] and |3] for the two high-density cells 
1500 Zr and 982 Ge. The proton density does not corre- 
spond to that of a proton cluster, as a non-negligible 
number of protons are sitting at the edge of the 1500 Zr 
cell, while the protons are spread out over the whole 
982 Ge cell. 

For the 1500 Zr cell, the proton density also depends 
on the strength of the pairing interaction (see bottom 
panel of Fig. [2]). For HF calculations (i.e., zero pairing 
strength) and for HFB calculations with Vj ow fc in the 
pairing channel, the box effect is much stronger than for 
calculations with the Gogny Dl and the DDDI pairing 
forces. 

These instabilities are sensitive to the Skyrmc func- 
tional used and to the boundary conditions as well. Neu- 
tron and proton HFB densities for the 1500 Zr cell are 
shown in Fig. [4] for calculations using the Vj ow k pairing 
interaction on top of a mean field built with SLy4, SkM* 
and MHF Skyrme functionals. The results obtained with 
the BCO boundary conditions are also shown for the 
SLy4 functional in the same figure. The proton instabil- 
ity in 1500 Zr is present only for SLy4 with BCE bound- 
ary conditions. All calculations for the higher- density 
cell 982 Ge give results very similar to each other, with 
the protons spread out over the entire volume of the WS 
cell. 

We conclude that solving the HFB equations in a WS 
cell of densities pb > 0.25po not always leads to a reliable 
solution. The method starts to be unstable at baryonic 
densities pb ~ 8 ■ 10 13 g • enr 3 « 0.25po an d it definitely 
breaks down at Pb ~ 10 14 g • cm -3 w 0.35po- 

We believe that the WS approximation should not be 
used whereas these instabilities occur. Hence, in the fol- 



lowing we show results only for the calculations where 
the HFB solution is stable. 



B. Pairing gaps 

It is useful to compare the superfluid properties of the 
inner crust to those of the INM. The pairing gaps at the 
Fermi surface are shown in Fig. [5] for INM and the inner 
crust. 

The lines represent the solutions A„(fc^„) of Eq. (|6]) 
in infinite neutron matter. The open points are the 
LCS pairing gaps obtained from HFB calculations in 
spherical boxes of radius Rws = 40 fm, with no pro- 
tons (Z=0) and with neutron number given by N — 
k Fn /(3n 2 ) -A-KR^g/Z. The solid points are the LCS neu- 
tron pairing gaps in the inner crust regions 11 through 2 
(see Table ©). 

As can be seen from the left panel of Fig. [51 the 
agreement between box and infinite-matter calculations 
is very good. However, the INM results obtained with 
the three pairing interactions differ substantially from 
each other. In particular, the gaps from the Vjowfc in- 
teractions are sensitively smaller than the others, with 
a maximum of 2.6 MeV against 3.2 MeV for the effective 
pairing interactions. At saturation, the Gogny pairing 
gap (~ 1.2 MeV) is much larger than the gap obtained 
with V\ owk (~ 0.2 MeV). The results with the DDDI 
pairing force differ even more. 

On the one hand, because of its simplicity, this type 
of contact pairing force is widely used in BCS and HFB 
calculations in INM and in the inner crust. Its parame- 
ters are usually fitted to reproduce given infinite-matter 
pairing gaps, and then used in inner-crust calculations. 
Hence, many parameter sets and energy cutoffs have been 
used in the past. On the other hand, realistic pairing 
interactions arc phasc-shift-cquivalent interactions, they 
require no fitting procedures and allow us to connect the 
theory to the basic nucleonic forces. The RG evolution to 
low-momentum and the separable representation in the 
channel make full HFB calculations in the inner crust 
feasible and fast, even on a laptop computer. We then ad- 
vocate for adopting these realistic pairing interactions to 
get more reliable results when applying microscopic the- 
ories to the inner crust of neutron stars. More than that, 
higher-order pairing correlations and other contributions 
to the pairing channel (i.e., three-body forces) represent 
necessary contributions [1, H3, HH and the adoption of 
realistic nucleonic potentials allows one to include them 
consistently. 

The effect of the proton clusters on the inner-crust 
pairing gaps is shown in the right panel of Fig. [5j The 
presence of the protons leads to a reduction of the pairing 
gap of at most 10% at the maximum of the curve. This 
effect is negligible below kp, n = 0.5 fm -1 . 

The pairing gap at the Fermi surface depends on the 
Skyrmc functional used, which defines the level density. 
Except for the highest-density regions close to the star 
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FIG. 1: HF (dotted lines) and HFB (solid lines) neutron and proton densities for the inner-crust regions 11 through 3 (see 
Table 10). The HFB results have been obtained using the Vi ow fe pairing interaction. The neutron (proton) density corresponds 
to the upper (lower) curves of each panel. 




FIG. 2: (Color online) Neutron (top panel) and proton (bot- FIG. 3: (Color online) Same as Fig. El for the 982 Ge WS cell, 
torn panel) densities for the 1500 Zr WS cell. The inset shows 
the proton densities in a semilogarithmic scale. 



6 




FIG. 4: (Color online) Neutron (top panel) and proton (bot- 
tom panel) densities for the 1500 Zr WS cell, using the VWfc 
pairing interaction on top of different Skyrme functionals and 
for different boundary conditions. The inset shows the proton 
densities in a semilogarithmic scale. 



core, there is a very small dependence on the bound- 
ary conditions and on the particular WS-cell (RwSiZ) 
paramctrizations. The pairing gap obtained with the 
Mow k pairing interaction on top of the Skyrme function- 
als SkM*, MHF and SLy4 is shown in the left panel 
of Fig. [6] The INM gap value for the three function- 
als reflects the different INM effective masses m*/m„ at 
saturation density (namely, 0.7531, 0.8687 and 0.997 for 
SLy4, MHF and SkM* respectively). 

The dependence of the pairing gap on the boundary 
conditions (namely, BCE and BCO) is shown in the right 
panel of Fig. HI The gap s obtained using the {Rws i Z) 
parameters from Ref. [20j are also shown in the same fig- 
ure. All points lie in a narrow band. One always has to 
keep in mind that either Refs. [3,[2(| used functionals and 
pairing interactions that differ from the Skyrme function- 
als and from the pairing interactions used in this work. 
One should use WS-cell (Rws,Z) parameters obtained 
from an energy-minimization procedure performed using 
the same functionals and pairing interactions later used 
to study the superfluid properties of the system. While 
this goes beyond the goal of the present work, we checked 
that even Ref. [2(| predicts a 10% suppression of the pair- 
ing gap at its maximum due to the presence of the pro- 
tons. 



C. Spatial extension of the Cooper-pair 

This section is dedicated to the study of the spatial 
properties of the Cooper pairs. The Cooper-pair wave 
function is defined as an [/^-weighted superposition of 
two-particle-state wave functions [0(ri)</>(r2)]oo coupled 
to total angular momentum J = and total spin S = 



$«(r 1 ,r 3 )= J2 



2j + l 



(ri)0(r 2 )] OO . (9) 



inn'lj 



In the following we will show the Cooper-pair wave 
function in the center-of-mass reference frame. There 
is a weak dependence of $ 9 on the angle between the 
Cooper-pair center of mass R and the relative position 
of the two nuclcons r. Hence, an angular average can be 
performed without loss of information. Results for the 
neutron wave function r) in the 1350 Sn WS cell are 

shown in Fig. [7] The dependence on the center of mass 
vanishes outside the nucleus, for R > 10 fm, recovering 
the asymptotic INM behaviour for large r values 

*W') « \k (r/ir Xp ) sin(A*, n r), (10) 

where Kq is the modified Bcssel function and \ r p is the 
neutron Pippard coherence length, defined as 



Xp 



h 2 kF,n 

m* n irA F ' 



(11) 



The modified Bessel function, whose expression for large 
values of r is K Q {r/-KX P ) ~ {x P /r) 1/2 exp[-(r/irXp)h 
damps the oscillating behavior of the INM coherence 
length. 

From the Cooper-pair wave function $™(i?, r), one can 
extract the coherence length. The coherence length gives 
information about the spatial extension of the pair. In 
INM, this quantity can be simply approximated by the 
Pippard formula defined in Eq. 1111 within an acceptable 
level of accuracy. See Ref. [36| for a detailed discussion. 
For a given density, the Cooper-pair extension in INM 
is uniform and its value is inversely proportional to the 
pairing gap at the Fermi surface. At variance, in finite 
nuclei, the neutron coherence length depends on the dis- 
tance R from the center of the nucleus and it is defined 



as 



x n {R) 



>fd 3 r r 4 |$"(.R,r)| 2 
Jd 3 r r 2 \<S> n (R,r)\ 2 ' 



(12) 



The neutron coherence length for Z = 40 and Z = 50 
systems is shown in Fig. [5] for a wide range of the neu- 
tron number. In the case of finite nuclei (i.e., 120 Zr and 
120 Sn), there is a minimum at the nuclear surface, as 
recent studies already discussed [37l - l40| . At the inner- 
crust densities this nuclear-surface effect gradually disap- 
pears as one moves from the outermost layers toward the 



7 




FIG. 5: (Color online) Pairing gaps obtained with the SLy4 Skyrme functional and the three pairing interactions in the legend. 
The lines represent the pairing gaps in infinite neutron matter and correspond to the solution A n {kF,n) of Eq. ([6]). Diamonds, 
squares and circles correspond to box calculations with DDDI, Gogny Dl and Vi ow fc pairing interactions, respectively. The 
open points in the left panel are the LCS pairing gaps in infinite neutron matter, obtained from spherical-box HFB calculations 
with no protons (see text). The solid dots in the right panel are the LCS neutron pairing gaps in the inner crust, for the WS 
cells 11 through 2 (see Table (0). 
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FIG. 6: (Color online) Dependence of the pairing gap on the Skyrme functional, on the boundary conditions and on the WS-cell 
(Rws,Z) parameters. All results are obtained with the Vi ow fc pairing interaction on top of different Skyrme functionals. The 
lines represent the pairing gaps in infinite neutron matter and correspond to the solution A n (fcF,n) of Eq. ((6j). The points 
represent the pairing gaps for the inner crust, obtained in box calculations. The solid line and the solid dots are the same 
as in Fig. [5] and are obtained with the SLy4 Skyrme functional and BCE boundary conditions. Left panel: HFB solutions 
obtained with SkM*, MHF and SLy4 Skyrme functionals and with BCE boundary conditions. Right panel: results are obtained 
with the SLy4 Skyrme functional. Solid and graded dots correspond to the two different BCE and BCO boundary conditions, 
respectively. The striped dots are the HFB solutions for the WS-cell (Rws, Z) parameters of Ref. [20I ]. 



star core (i.e., from 180 Zr and 500 Zr). In the outer neu- 
tron gas (i.e., large R), WS-cell calculations recover the 
INM coherence length value at the corresponding den- 
sity, as it can be seen from the arrows in Fig. [5] These 
arrows correspond to the Pippard coherence length Xp 
computed with the Fermi momentum kp,n of the outer 
gas of the given inner-crust region (see Table (|l])) and 
with A„(fci? >n ) taken from the INM results in the left 



panel of Fig. [5] 

At present, the relation between the coherence length 
and the strength of the pairing correlations still needs to 
be better clarified. According to recent studies in finite 
nuclei [37H40| , the coherence length has a very small de- 
pendence on the strength of the pairing interaction and, 
consequently, on the value of the pairing gap at the Fermi 
level. This lead to the conclusion 39] that the minimum 
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FIG. 7: (Color online) Cooper-pair wave function for 1350 Sn 
as a function of the relative distance r of the two nucleons of 
the pair. Different curves correspond to different center-of- 
mass values R. 
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FIG. 8: (Color online) Neutron coherence length for different 
Z = 40 (top panel) and Z = 50 (bottom panel) systems, 
calculated using the SLy4 functional and the Viawk pairing 
interaction. The arrows on the right correspond to the INM 
coherence length taken at the Fermi momentum of the outer 
neutron gas (see text for details). 



of the coherence length in finite nuclei has little to do 
with an enhanced strength of the pairing correlations at 
the nuclear surface. The minimum is a finite-size effect. 

A possible explanation is that the superfluid properties 
are mostly determined by the levels around the Fermi en- 
ergy (typically in the region Sp ± A). At variance with 
the INM case, where a continuum of states in the region 
Ef ± A contributes to the pairing correlations, in finite 



nuclei close to the stability valley the dependence of the 
coherence length on the pairing gap at the Fermi level is 
washed out by the shell structure. The situation is differ- 
ent for the dineutron halo nucleus 11 Li [4l|. The weakly- 
bound valence single-particle wave functions extend far 
outside the core and couple with the continuum. As a 
consequence, finite-size effects are not able to suppress 
the dependence of the coherence length on the pairing- 
interaction strength. The coherence length in 11 Li has a 
minimum at the nuclear surface, but this minimum dis- 
appears with a negligible pairing interaction. 

In the same way as 11 Li, we can then expect the coher- 
ence length in the inner crust to depend on the pairing- 
interaction strength, as the Fermi energy lies in the con- 
tinuum. Following the ideas of [39|, lU |42| in finite nuclei, 
we investigated this dependence in the 180 Zr and 1800 Sn 
WS cells, whose neutron coherence length is shown in 
Fig.[5]for the SLy4 Skyrme functional and with the V\ OV! k 
pairing interaction rescaled by a factor x. The curves in 
each panel of Fig. correspond to different rescaling fac- 
tors x. From top to bottom, pairing correlations increase, 
with consequently larger and larger LCS pairing gaps. 
The coherence length in the inner crust depends on the 
strength of the pairing interaction. This dependence is 
present also inside the proton cluster and it is stronger in 
the outer neutron gas, where the inverse proportionality 
to the pairing gap is recovered (see Eq. fTTj) . 

We conclude that the coherence length can depend on 
the strength of the pairing interaction, even inside the 
nucleus, but finite-size effects suppress this dependence 
when the Fermi level docs not lie close to the continuum. 



D. Pairing field 

In this section we investigate the spatial properties of 
the pairing field. For a local pairing interaction u(ri— v^), 
the pairing field reads [13] 



A»(n,r 2 ) = -u(n -r 2 )$ 9 (r 1 ,r 2 ), 



(13) 



where $ 9 (ri, r 2 ) is the Cooper pair wave function defined 
in Eq. [9l It is convenient to perform a Wigner transform 
[2f| , changing to the center-of-mass coordinates and per- 
forming a Fourier transform on the relative distance r. 
In this way we obtain the pairing field A 9 (R, ki 2 ), as 
a function of the two-particle center of mass R and of 
their relative momentum ki2- The pairing field depends 
weakly on the relative angle between the vectors ki 2 and 
R, so we show the results after performing an angular 
average. 

The neutron pairing field A n (R, &i 2 ) is shown in 
Fig. Ql for the 1350 Sn WS cell, using the SLy4 functional 
and the Vj ow fc pairing interaction. We observe that the 
pairing field A n (R, fci 2 ), has a strong dependence on the 
relative momentum k\2 and it is suppressed at the center 
of the cell by the presence of the nucleus. As a conse- 
quence one observes a global reduction of the pairing gap 
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FIG. 10: (Color online) Neutron pairing field A n (R,k 12 ) for 
1350 Sn (solid lines) as a function of R for fixed values of the 
relative momentum fci2. The dashed line corresponds to the 
local neutron pairing field Al OC (_R) defined in Eq. 1151 The 
dot-dashed line corresponds to the LDA approximation [il ]. 



field agree 



(cf. Sect. ITnBj) . Our results on the pairing 
with the more detailed discussion of Ref. [43[ . 

A local approximation lA\ oc to the pairing field can 
be obtained by defining a local Fermi momentum kp^ q (R) 



2m q (R) 

and then taking the pairing field value 

A q (R, k F , q {R)) 



(14) 



(15) 



Ucentr(R) 1S the central potential of the single-particle 
Hamiltonian. 

The neutron local pairing field A£ oc is shown in 
Fig. [TU1 together with an LDA approximation Alda to 
the pairing field (44J. Alda has been obtained solv- 
ing the neutron gap equation (cf. Eq. [5]) in asymmet- 
ric nuclear matter, at different values of the Fermi mo- 
mentum kF, n (R) — (3vT 2 p„(i?)) 1 / 3 , where p n (R) is the 



HFB density for the 1350 Sn cell. The LDA pairing field 
is more suppressed in the interior of the nucleus (where 
kp.n ~ 1-33 fm -1 ). and it is peaked at the nuclear sur- 
face. The local pairing field, instead, has a monotonic 
behavior. The reason of the difference between the two 
approximations has already been analyzed in a previous 
work [45[ and comes from the fact that the classical LDA 
approximation is not able to capture the shell-structure 
effects of the system. 

The evolution of the local neutron pairing field from 
finite nuclei to the inner crust is shown in Fig. [TT] for 
Z = 50 systems. The analysis is done for both V| ow & (left 
panel) and the DDDI (right panel) pairing interactions. 
The latter is a local interaction and the pairing field docs 
not depend on ki2 by construction. As we go from finite 
nuclei to the inner crust, the role of the nuclear surface 
changes drastically. The local pairing field is peaked at 
the surface of finite nuclei, at variance with nuclei im- 
mersed in a sea of supcrfluid neutrons. The value of the 
pairing field inside the nucleus represents the main differ- 
ence between the results obtained with the two pairing 
interactions of Fig. [TT] 

Although the pairing field and its local approximation 
can be interpreted as a measure of the strength of the 
pairing correlations, one has to keep in mind that these 
are not observables. Hence, the association of a small co- 
herence length (i.e., close correlated nucleons, see Fig. [BJ) 
with a large pairing field (i.e., large pairing correlations, 
see Fig. [TT)) has to be taken with a grain of salt. The sup- 
pression of the dependence of the coherence length on the 
strength of the pairing interaction, due to shell-structure 
effects in finite nuclei makes the above association at least 
unclear. 



IV. CONCLUSIONS 

We studied the supcrfluid properties of the inner crust 
of neutron stars, representing the different density re- 
gions with 11 spherical WS cells. In the innermost layers 
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FIG. 11: (Color online) Neutron local pairing field for Z = 50 systems, obtained with VWfc (left panel) and the DDDI (right 
panel) pairing interactions on top of the SLy4 Skyrme functional. 



close to the star core, the WS approximation turned out 
to break down, even at the HF level, where protons leak 
out of the center of the cell. The method starts to be 
unstable at baryonic densities pb ~ 8 • 10 13 g • cm" 3 and 
it definitely breaks down at pb « 10 14 g • cm -3 . The WS 
approximation should not be used whereas these insta- 
bilities occur. 

Within the limits of applicability of the method, we 
performed fully self-consistent HFB calculations, based 
on Skyrme functionals plus realistic pairing interactions. 
The pairing interaction has been taken as the low- 
momentum evolution of the Argonne AV18 potential. We 
restricted our calculations to the 1 So pairing matrix el- 
ements, as higher partial waves are expected to give a 
much smaller contribution in the inner crust. 

From a comparison with INM, the presence of the pro- 
tons in the inner crust is found to reduce the pairing gap 
at the Fermi surface of about 10% at its maximum. This 
suppression is negligible below kp^ n = 0.5 fm _ . We also 
compared two different WS (Rws, Z) configurations, ob- 
tained from two energy-minimization procedures [H, l20| , 
with the results lying in a very narrow band. 

The results obtained with effective pairing interac- 
tions (namely, the Gogny Dl interactions and a density- 
dependent contact force) differ substantially from the re- 
sults obtained with the realistic pairing interactions. Not 
only the non-arbitrariness of the latter ensures a higher 
reliability, but also the availability of a high-precision 
separable representation of low-momentum realistic po- 
tentials makes the calculations feasible even in the inner 
crust of neutron stars. Hence, we advocate for adopting 
these realistic pairing interactions to get more reliable 
results when applying microscopic EDF theories to the 
inner crust of neutron stars. 

Higher-order pairing correlations are expected to play 
an important role in the inner crust of neutron stars, 
where the exchange of collective vibrations leads to a re- 
pulsive interaction and, consequently, to a suppression 
of the pairing gap [HI, S] . At sufficiently high densities 



(&F,n > 0.7 fm -1 ), three-body forces are also expected 
to contribute to the pairing interaction 47J. A quanti- 
tative assessment of this effect requires all contributions 
to be treated consistently and represents the subject of a 
future work. Recently, ab-initio Quantum Monte Carlo 
calculations have been carried out for low-density infinite 
neutron matter (kp^ n < 0.5fm _1 ) (48|, with an effort to 
understand the differences with other Monte Carlo re- 
sults (49l - [5l| . A quenching of the mean- field pairing gap 
is predicted, resulting in a larger pairing gap than what 
previous calculations found [52h56| . 

In this work we also presented a detailed study of the 
Cooper-pair spatial properties, which showed how the co- 
herence length in the inner crust depends on the strength 
of the pairing interaction. This dependence is present 
also inside the proton cluster and it is stronger in the 
outer neutron gas, where the inverse proportionality to 
the pairing gap is recovered (see Eq. fTTj) . This result was 
expected at the inner-crust densities, where the Fermi en- 
ergy lies in the continuum, at variance with nuclei close 
to the stability valley, where shell-structure effects sup- 
press this dependence. At present, the relation between 
the coherence length and the strength of the pairing cor- 
relations still needs to be better clarified. 
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